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It is shown how any Lindbladian evolution with selfadjoint Lindblad operators, either Markovian 
or nonMarkovian, can be understood as an averaged random unitary evolution. Both mathematical 
and physical consequences are analyzed. First a simple and fast method to solve this kind of master 
equations is suggested and particularly illustrated with the phase-damped master equation for the 
multiphoton resonant Jaynes-Cummings model in the rotating-wave approximation. A generaliza- 
, tion to some intrinsic decoherence models present in the literature is included. Under the same 

' philosophy a proposal to generalize the Jaynes-Cummings model is suggested whose predictions 

' are in accordance with experimental results in cavity QED and in ion traps. A comparison with 

, stochastic dynamical collapse models is also included. 
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I. INTRODUCTION 

•/^ ■ Since the early years of quantum mechanics W the principle of quantum superposition has been recognized to play a 
prominent role in the theory and its applications. The destruction and preservation of these superpositions of quantum 
states occupy a central place in issues such as the quantum-to-classical transition ^ and potential technological 
applications in Quantum Information, Computation and Cryptography |[ ^. From a physical standpoint the 
, loss of coherence in quantum systems is rooted on the pervasive action of the environment upon the system. This 

■ environmental action has received a careful mathematical treatment (cf. |0, ^, |^ and multiple references therein) 
I going from a constructive approach based on disregarding the degrees of freedom of the environment due to their 

lack of control by the experimenter (" tracing-out" methods) to an axiomatic approach based on the initial setting of 
physically motivated axioms to derive an appropiate evolution (master) equation for the system [l^ . 
, Most of these master equations (ME's hereafter) satisfy the Markov approximation (semigroup condition) and can be 

■ put into the Lindblad form : 

^ ; 
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^ = -^[H,p{t)] + ^Y.[[V,p{t),VJ] + [V„p{t)vj} (LI) 

j 

^ . where H is the Hamiltonian of the system and {Vj} are operators (so-called Lindblad operators) containing the effect 
■ - ' of the environment upon the system. Indeed in the axiomatic approach the Markov approximation is posed as an 
initial hypothesis [ pO[ , thus rendering highly difficult a generalization to nonMarkovian situations. 

In this work we develop a novel attempt to derive ME's both in the markovian and the nonMarkovian regimes using 
stochastic methods ^ jointly with well-known operator techniques commonly used in quantum mechanics [Q. 
The main idea consists of building random evolution operators (evolution operators with one or several stochastic 
parameters in it) which contains the decohering effect of the environment and then taking the stochastic expectaction 
value with respect to this (uncontrollable) randomness. The paper is organized as follows. In section || we state 
and prove our main (though still somewhat partial) result, namely that any Lindblad-type ME, either Markovian or 
nonMarkovian, with selfadjoint Lindblad operators can be understood as an averaged random unitary evolution. In 



section [II we discuss some first mathematical consequences of this result such as a very fast method to solve ME's 



provided the unitary solutions are known; we illustrate this by solving the phase-damping ME for the multiphoton 
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resona nt Jay nes-Cummings model in the rotating-wave approximation 
section 



a (section [II A). We then comment in 



III B two immediate consequences, namely both Markovian and nonMarkovian regimes are attainable under 



the same mathematical formalism and the Lindbladian structur e with selfadjoint Lindblad operators is shown to have 
an origin independent of the Markov approximation. In section III C we show how the flexibility of the mathematical 
language employed can easily generalize some intrinsic decoherence models present in the literature In 
section |^ we discuss the previous ma in res ult under a more physical spirit by proposing a slight generalization of 
the Jaynes-Cumming s mo del (section IV A), c ompa ring this proposal with experimental results in optical cavities 
experiments (section [V B ) and finally (section IV C ) showing how the proposed formalim can account for reported 
exponential decays of Rabi oscillations in ion traps. We include in section ^ some important comments regarding a 
brief comparison with existing models of stochastic evolution in Hilbert space, the possibility of intrinsic decoherence 
phenomena and future prospects. Conclusions and a short appendix close the paper. 



II. LINDBLAD EVOLUTION AS AN AVERAGED RANDOM UNITARY EVOLUTION 

The main result whose consequences are to be discussed below is the following: Every Lindblad evolution with 
selfadjoint Lindblad operators can be understood as an averaged random unitary evolution. We will analyse this 
proposition in detail. The objective is to reproduce the Lindblad equation 



dpjt) 
dt 



n 

= -i[H, Pit)] + - J2mpit),v^ + m,pm]} 

i=l 
1 " 

= -^[H,p{t)]--J2[V^,[V,,pm 



(11.2) 



by adequately modifying chosen parameters in the original evolution operator. For simplicity let us start by considering 
the case n = 1. We will first study the case where the Hamiltonian H and the (selfadjoint) Lindblad operator V 
commute. It is very convenient to introduce the following notation. The commutator between an operator G and X 
will be denoted by Cci-'^] = [G, X]. Thus the von Neu mann -Liouville operator will be £ = —iCn, where H denotes 
the Hamiltonian (h — 1). Then the Lindblad equation (11.2) with n = 1 can be arrived at by 



1. Adding a stochastic term BtV to the argument of the evolution operator: 



U{t) = exp{-itH) ^ Ustit) = exp{-itH - iBtV) 

where Bt denotes standard real Brownian motion [ p^ . 
2. Taking the stochastic average with respect to Bt in the density operator deduced from Ust{t)'- 



(II.3) 



p{t) = nust{t)p{mlm 

where E denotes the expectation value with respect to the probability measure of Bt- 



(II.4) 



The proof of this result is nearly immediate. Taking advantage of the commutativity of H and V and making use 
of theorem 3 in (cf. appendix; relation (A.3)) we may write for the density operator: 



p{t) = exp(i/:)E[exp(-i6tCy)][p(0)] (II.5) 

Thus all we have to do is to calculate the expectation value of the random superoperator exp(— iStCy). Developing 
the exponential into a power series and recalling = -pr^i" if is even and IE[.B"] — otherwise, we arrive 

at 



pit) = exp(-f/:)exp(-^C2,)[p(0)] 
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exp(. 



-tC--Cl)[pm 



(11.6) 



which produces the desired master equation: 



dp{t) 
dt 



^i[H,p{t)]~-[vAv,pm 



(II.7) 



When the Hamiltonian H and the Lindblad operator V do not commute, the previous method is not suitable, since 
the calculation of the expectation value cannot be performed in the same way. A way to circumvent this problem 
is to proceed in the same way as before but in the Heisenberg picture. Thus let U{t) = exp{—itH) be the original 
evolution operator in the Schrodinger picture. The corresponding evolution operator in the Heisenberg picture will 
trivially be Unit) — I- As before we proceed in steps: 

1. Add a stochastic term to the argument of the evolution operator: 

Unit) = eM-itO) ^ UHM(t) ^ 

= Texp(-z / VH{s)dB,) (II.8) 
Jo 

where T denotes time-ordering and Vnit) is the operator V in the Heisenberg representation. 

2. Take the stochastic average with respect to Bt in the corresponding density operator: 



PHit)^nuH,stit)pm}ist] 



(H.9) 



3. Finally to arrive at (H.2) change to the Schrodinger representation. 



A comment should be made. T he st ochastic term added to the original evolution operator in (H.8) is a natural 
generalization of the one added in (H.3). The Ito integration now appears as a consequence of the time dependence 
of the operator to be added: (i). Note that when [H, V] =0, Vff(t) = V and the stochastic term reduces to VBt 
as before. 



Now to perform the previous tasks is a bit more involved. Firstly combining relation ( A. 3 ) and the T— operation, the 
step |l| can be carried over: 



(H.IO) 



The expectation value (H.IO) can be evaluated by resorting to functional techniques [T^. Recall that the 
characteristic functional of a stochastic process Xt is defined as 



exp i I k{t)xtdt 



where k{t) is an arbitrary real- valued function. In particular, for white noise Xt — £,t (cf. |lq|) 



(11.11) 



G^[k{t)] ^ E[exp i / k{s)dB 



where dBt = Ctdt, E^t = and E[^t^s] = 5{t — s) have been used. From this it is then clear that ( 11.10 ) can be written 
as 



(11.12) 



PHit)=Te---^oCl„is,ds^p^Q^^ 



(H.13) 
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Back to the Schrodinger picture, the master equation derived from (11.13) is 



dpjt) 
dt 



-i[H,p{t)]--[V,[V,p{t)]] 



(11.14) 



The generahzation to many Lindblad operators is elementary: aU we have to do is to use the n-dimensional standard 
real Brownian motion ||l^ Bt — {B}, ■ ■ ■ ,B^). The strategy is the same. 

III. ANALYTICAL CONSEQUENCES 

The first consequences one can derive from the previous result are of analytical fashion. As an immediate aplication 
we will show how the Jaynes-Cummings model with phase damping in the rotating-wave approximation can be solved 
provided we know the solution to the original Jaynes-Cummings model. As a second consequence we will discuss 
how the previous result can be generalized to nonMarkovian situations, thus providing a common language for both 
Markovian and nonMarkovian evolutions. Finally it is shown how existing intrinsic decoherence models are naturally 
generalized using this formalism. 

A. Solution of the Resonant Multiphoton Jaynes-Cummings Model with Phase Damping 

The Jaynes-Cummings model (JCM hereafter) |2l] shows an undoubtable relevance in the study of quantum 
systems in different fields such as Quantum Optics, Nuclear Magnetic Resonance or Particle Physics. It is an exactly 
solvable model which allows us to study specifically quantum properties of Nature such as electromagnetic field 
quantization or periodic collapses and revivals in atomic population. The JCM describes the evolution of a two- 
level quantum system (the atom) interacting with a mode of the electromagnetic field under certain approximations 
(rotating wave approximation, dipole approximation, etc. - cf. ||2^, |2l] for details). Usually in normal experimental 
conditions this will be an idealization and the environment should be taken into account, the effect of which can be 
very appropiately treated introducing a phase-damping term Q . Thus the master equation for this system will read 



^ = -^[H,p{t)]-l[H, [H,p{t)]] 



(III.15) 



where H is the Hamiltonian of the system and 7 is a damping constant. Here we will show how ( III.15| ) can be very 
easily solved when H is the resonant multiphoton JC Hamiltonian (cf. ||2^ for an alternative approach), i.e. when 



H = ija^a + LJoS, + A(S'_a^™ + S+a") 



(HI. 16) 



where u denotes the frequency of the field mode, ujq is the atomic transition frequency, A is the atom-field coupling 
constant, a' and a are the mode creation and annihilation operators respectively, Sz is the atomic-inversion operator 
and S± are the atomic "raising" and "lowering" operators. An exact resonance is assumed, thus loq = muj. 
We will focus in two quantities of relevant physical meaning, namely the atomic inversion W{t) = TT[p{t)Sz] and 
the photon number distribution at time t: Pn{t) = {n\TrAp{t)\'n) . To compare with methods found in the literature 
p2| , p3| we will restrict to the case in which initially the atom is in its excited state |-|-) and the electromagnetic field 
is in a coherent state \a) = Yl^=aQnW) ^ with Qn = exp(— |ap/2)^^. The unitary evolution (7 = 0) provides the 
following expressions for these quantities: 



Tl = 

Pn{t) = IQnpCOS^ \t 



2Xt 



|Q„_mpsin^ 



Xt 



(n + to)! 



{n + my. 



[n + to)! 



(IH.17a) 



(IH.17b) 
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The objective is to calculate these same quantities when the phase damping term is present in ( III.15|), i. e. when 
7^0. We will make use of the result proved in the previous section and note that the equation ( [III. 15 ) can be 
obtained by adding a stochastic term to the original evolution operator and then performing the stochastic average. 
In our case, V = , which obviously commutes with the Hamiltonian, thus we are in the first case. The original 

evolution operator is promoted to 



U{t) = exp{-itH) 



Ust{t) = eM-itH - i-f^'^BtH) 



(III.18) 



Equivalently we may think that it is t which is promoted t ^ t + j^/^Bf . Thus to arrive at the desired "phase- 
damped" expressions WP'^{t) and Pn'^{t) all we have to do is to add a stochastic term to the time variable t and then 
perform the average: 



,n=0 



(III. 19a) 



IQriPcOS^ 



|(3„_™psin^ 



(III.19b) 



Using the linearity property of the expectation value and recalling the moments of the standard real brownian 
motion (cf. above and [|l2|), the previous calculations can be carried over elementarily using (see appendix ^): 



Hence 



E cos 



[2A(, + ../»R,\/^l 



X COS 



2Xt 



(n + to)! 



(III.20) 



2,(n + m)! 



COS 



2Xt 



(n + m)! 



n! 



1 + e 



I COS 2\t\j 



{n + my. 



^|Q„-™P{l-e-2-^^^*^}cos 



2Xt 



{n + my. 



(III.21a) 



(III.21b) 



which exactly coincides with equations (41) and (43) in lE2| and eqs. (3.26) and (3.27) in [gsj for to = 1. We encourage 
the reader to compare this method with those used in 12^ |2^ . 

Obviously this formalism can also be used to solve the equation ( III.15| ) with any arbitrary Hamiltonian provided we 
already know the solution when 7 = 0. 
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B. NonMarkovian Evolution 



A second consequence of the formalisni depicted above is its immediate generalization to nonMarkovian situations. 
The resuh in section || can be readily generalized to the following: Any Lindbladian master equation, whether 
Markovian or nonMarkovian, but with selfadjoint Lindblad operators can be obtained as the stochastic average of a 
random unitary evolution. The generalization stems out from the single fact that whereas in the Markovian regime we 
necessarily have to add a stochastic term of the form BtV — V /p dBt, in the nonMarkovian case this restriction drops 

out and then we may add a term like V v{s)dBs, where v{s) is an arbitrary real-valued function which encodes e.g. 
the time response of the environment to the system evolution p^ . Under these circumstances, the previous procedure 
(for simplicity's sake we will only care about the commuting case; the noncommuting case is similar) drives us to 



p(t) = exp(t£)E[exp ^-iCv^ vis)dBs^][p{0)] (III.22) 

Now developing the exponential again into a power series, calculating the expectation value of each term with some 
elementary Ito calculus and resumming the series, one arrives at 



^ ^ -^[H,p{t)] - 2^[V, [V,pm (III.23) 

where ^{t) — w^(s)(is. This is clearly a Lindbladian nonMarkovian master equation. The extension to more than 
one Lindblad operator is again trivial. This result casts some light into the origin of the Lindbladian structure of 
master equations with selfadjoint Lindblad operators, independently of their Markovian or nonMarkovian character, 
something beyond reach of the original axiomatic approach of [|lO[ |ll| . 

In this sense the result proven here generalizes previous derivation of Lindblad evolution using stochastic calculus 
Ip^ , |25|| by dropping out the semigroup condition. Note that this generalization allows us to conclude that since 
7(t) = v'^{s)ds the decoherence process is irreversible, i.e. no coherence can be recovered within the domain of 
validity of the phase-damping ME as an evolution equation for the quantum system. 

The time dependence of the decoherence factor also suggests a classification of different kinds of environments de- 
pending on the rate at which the environment decoheres the system (cf. ^^). It remains open what the physical 
conditions should be to have the different decoherence factors. 



C. Models of intrinsic decoherence 



A third advantage appears as a natural generalization of intrinsic decoherence models already present in the liter- 
ature p6| . These two models propose an intrinsic mechanism of decoherence based on the random nature of time 
evolution (we will not enter into the discussion of the physical justification of this hypothesis -see original references 
f or disc ussion, we will only show how they can be generalized), which basically drives us to the evolution equation 
( 111.15 ). The starting hypotheses (apart from the random nature of time evolution and the usual representation of a 
quantum system by a density operator) are a specific probability distribution [l^ or a semigroup condition (Marko- 
vianity) |17| for the time evolution. As a result we obtain a nondissipative Markovian master equation in both cases. 
The formalism presented here dispenses with any of these specific conditions, something which allows us to obtain 
more general master equations, i.e. both Markovian or nonMarkovian and dissipative or nondissipative. 
The result comes from the combination of Ito calculus and the spectral representation theorem for unitary operators 
|p7| . Let U{t) be an unitary evolution operator. By means of the spectral decomposition theorem it can be 
written as 



U{t) = I e-'^'dEg (III.24) 



where Eg denotes the spectral measure of the evolution operator. Now we perform the stochastic promotion as before 
by substituting 

Ot~^XtiO) (in.25) 
where XtiS) is real stochastic process. Then (see 12^ for details) 
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1. The Markovian nondissipative master equation appearing in Q and Q is obtained if Xt(^) = Ot + j^/'^Bt- But 
note that now the Lindblad operators are not fixed by any initial assumption. If e.g. Xt(^) = Jq ^{^'^ d)dBs{0) 
with correlation function dBt{0)dBt{9') = e-'^^i^-^ f dt we arrive at a Lindblad equation 



dt 

2 



(III.26) 



which is clearly different from the phase-damping master equation (III. 15). Thus even restricting ourselves to 



the same range of assumptions (Markovianity and nondissipation) we can obtain more master equations. 

2. A nonMarkovian nondissipative master equation like e.g. 

^ = -^[H,p{t)] ^[H, [H,pit)]] (III.27) 

is obtained if XtiS) = Ot + /J a{s)dBs and \{t) = /J <7^{s)ds. 

3. More general equations can be readily obtained by appropiately combining the correlation properties of Bt{0) 
and the time dependency of a{t; 9). 

This allows these models to be used to explain a wider range of phenomena than that originally considered. 



IV. RABI OSCILLATIONS DECAY IN CAVITY QED AND ION TRAPS EXPERIMENTS 

In previous sections we have developed a method to adequately modify the original evolution operator of a quantum 
system to finally arrive at a Lindbladian master equation. Now we find legitimate to proceed the other way around, i.e. 
what physical predictions are derived from the assumption that a parameter in the evolution operator of a quantum 
system is random? To be concrete we will focus upon two different physical systems, namely a Rydberg atom in an 
optical cavity and a linear rf (Paul) ion trap. We will confront the previous hypothesis with experimental results. 



A. The Jaynes-Cummings Model Revisited 



The JCM model describes the interaction between an atom and the electromagnetic field under very special condi- 
tions |2^, Different generalizations have been proposed to take the model closer to experimental reality keeping 
its solvability. Among these one can find the inclusion of dissipation (often modelled by coupling the field oscillator 
to a reservoir of external modes) and/or damping (as a consequence of spontaneous emission), multi-atom, multi-level 
atom, generalized- interaction and multiple-mode generalizations (see |2^ for references). 

Here we want to introduce a novel proposal, which states that JCM predictions can be rendered more realistic by 
noticing that the coupling constant A between the atom and the field mode should have a stochastic part which 
contains part of the effects of the approximations assumed in constructing the original model. Since these effects are 
not under control, to make physical predictions we must average on the introduced random parameters. To illustrate 
the idea let us consider the original JCM with Hamiltonian H = Loa^a -\- ujqSz + \{S-a) -\- S'+a), where lo denotes 
the frequency of the field mode, a and a) their corresponding creation and destruction operators, wq the frequency 
difference between the two energy atomic levels, Sz the atomic population operator, A the atom-field coupling constant 
and S± the energy raising/lowering atomic operators. We claim that the evolution stemming out from H should be 
modified by inserting a random part Xt? where Xt is a real stochastic process which contains the departure from the 
original ideal situation. The connection with the previous formalism is established by noticing that the evolution 
operator (in interaction picture) must then be: 



Ui{t) = T exp 



{\ + Xs)HT\s)ds 
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exp 



^i{Xt+ / Xsds)H'' 
Jo 



(IV.28) 



where = {S-a^ + S+a) is the interaction Hamiltonian and H}^^{t) = ul{t)H"^^Uo{t) is the interaction Hamil- 
tonian in interaction picture (for simpHcity's shake exact resonance has been assumed oj = ujq). Now the expression 
Af = At + /q Xsds is a real stochastic process itself which can always be expressed in the form 



At = E[At] + / v,dB^ 



(IV.29) 



where Vt is a real stochastic process uniquely determined by Af. Then the density operator in interaction picture will 
then be given by 



Pi{t) 



cxp 



(E[A,] + VsdB,) Ch 



ds 



(IV.30) 



Instead of giving the general form of the expectation value in ( IV.3C| ) (which will be difficult to obtain in full 
generality), we will propose some physical choices based on xt- If Xt = I'^^'^^t^ i-e. the original deterministic 
evolution is randomly perturbed by a white noise coming from a stochastic perturbation of the coupling constant, 
then \t^\t + 'f^^^Bt and (|IV.30D reduces to 



pi{t) = exp{~tXtCH,)E[cxp(^-ij^/''BtCH,)[p{Q)] 
= E[cxp(-*(Ai + 7i/2St)CH,)]b(0)] 

which yields a density operator in Schrodinger picture given by 



(IV.31) 



pit) = E[exp (^-itCHo -i{xt + -i^l^Bt) Ch,)][p(0)] (IV.32) 

This relationship means that to obtain the physical predictions of this proposal, all one has to do is to make the 
sustitution At ^ At + j^^^Bt in the original expressions and calculate the expectation value. 

Note that this proposal allows us to embrace nonMarkovian (though Lindbladian) situations with little extra effort, 
e.g. by claiming that the random perturbation is time-dependent Xt = l^^'^{t)£,t- More general options are also 
possible. The particular choice for xt relies upon the specific system under study. Notice also that the generalization 
proposed here is compatible with the ones quoted above, i.e. one may combine both type of generalizations. 



B. Decay in an Optical Cavity 



Let us consider the situation depicted in [2^, which appears as the first direct (in time domain) experimental 
evidence of field quantization. The system consists of a Rydberg atom in a high-Q optical cavity with the atom 
initially excited and the electromagnetic field in a coherent state |a). This system is accurately described using the 
JCM, thus Rabi oscillations are expected and concordantly experimentally measured. The theoretical prediction for 
the probability Peg{t) to find the atom at a time t in its ground state is Peg{t) = IQ^P [Xty/n + l) . However 

an exponential damping in these oscillations are detected (see for details). The stochastic JCM accounts for this 
damping (see fig. 0) assuming xt = 7^^^Ct which produces 



P^g'^'^it) = ^ 1 - E IQ«pe-2("+i)^^* cos (2A^^A^TT) (IV.33) 

V n=0 I 

The physical interpretation under this assumption is rather clear: the ideal coupling assumed in the original JCM 
does not hold any longer and departures from this ideality should be considered. Darks counts and decoherence 
caused by collisions with background gas have been considered as candidates to explain the damping behaviour |45| 
and even more radical proposals appear in the literature . Except for the latter which reveals an original intrinsic 
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process, all of them resort to external agents. Here note that we do not need to do so, since the departure from ideal 
atom-field mode coupling can be justified within the domain of the JCM assumptions themselves, i.e. the assumption 
of coupling between a unique field mode and two levels of the atom can be relaxed by adding in a natural way a 
random background in this coupling. 




FIG. 1: Peg{t) for the original and the stochastic JCM's with parameters A — 50TvkHz, 7 — 1/2-k and |a|^ = Uphotons ~ 0.4. 
See H. 



This decay is also present in the case of arbitrary initial conditions, 
orthodox prediction, then the previous recipe drives us to 



If Pegit) = Er=o^"Sin2 {XtV^TTT) is the 



psJCM 



it) 



-2(n+l)7At 



cos 



{2XtV 




(IV.34) 



where P„ depends on the actual initial conditions of both the atom and field mode. Note that not only can any actual 
exponential decay (with arbitrary time dependency) be obtained with the substitution Xt Xt+ j^^^{s)dBs, but also 

possible changes in the argument of the cosine function could be accounted for by making Xt uj{t) + ^^l'^{s)dBs- 
See appendix ^ for details. 

In this way we have obtained a decohering system (oscillations coming from quantum superpositions are progressively 
supressed) without necessarily resorting to the action of the environment and keeping quantum principles untouched 
(see discussion in section ^ later on) . This opens new possibilities to discuss possible sources of decoherence. 

C. Decay in an Ion Trap 

The previous experimental supression of quantum coherence has also been detected in a linear Paul ion trap 
The physical situation is formally similar to that of the Rydberg atom coupled to a field mode: the laser 



31 



field is operated upon the trap in such a way that it can be considered that only two internal energy levels of the 
ions are coupled to the center-of-mass (COM) mode of the set of ions (see [0). In the dipole and rotating- wave 
approximations, the interaction Hamiltonian (in interaction picture) is 



iJ}"*(t) = f^5+exp [i [q{ae' 



-Ka^e*'"^*)-(5i]) 



(IV.35) 



where 77 = fczo is the Lamb-Dicke parameter ik = 2'K/Xiaser and zq = ((0|z|0)^/^), S± denote the raising/lowering 
operators for the internal levels, a and a) denote the destruction and creation operators for the COM mode, u)z 
denotes the frequency of the harmonic trap for the COM, 5 = lo — ojq with uj the frequency of the laser mode and ojq 
denotes the difference between the two internal energy levels of the ions. 
The statevector can then be written as 



® \n) 



(IV.36) 
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where \mz) and \n) denote the (time-independent) internal and motional eigenstates. In the conditions of interest, 
i.e. in resonant transitions {5 = LOz{'n' — n) with n,rt' integers), the coefficients Cm^,n{t) satisfy the equations |^l| 



(IV.37) 
(IV.38) 



where is given by ^n,n' = (n'le'*^^ |n) | (fc = 27r/Aiaser and z = zo(a + a^^) is the operator for the COM motion). 

From these one can predict the well-known Rabi oscillations of the system. For concreteness' shake let us focus upon 
the first blue sideband case, i.e. when n' = n + 1. If the trap is prepared in the initial state |5'(0)) = |— ) |«) = \ — 
then the probability of finding a single ion in the |— ) state at time t is P_ (i) — cos^(f2„,„+ii). However experimentally 
an exponential decay is obtained. As before, one can argue that ideality should be restricted and both a substitution 
^n.n+it — > fln,n+it + Jq ^^^'^{s)dBs and the corresponding averaging should be performed on P-{t). This would drives 
us to the relation 



P 



■rand 



it) 



-2A(t) 



cos(2fi„,„+ii) 



(IV.39) 



where A(t) = 'y{s)ds. This way the exponential decay would have been obtained. Physically this recipe can be 
justified by taking into account intensity fluctuations in the laser field (see -notice that some generality is gained 
with respect to this work). 

However experimental data for the COM initially in an arbitrary state and the ion in the ground state |— ) are better 
fit by 



cos(217„,„+it) 



(IV.40) 



where P„ is an n-dependent quantity which relies upon the initial conditions of the ion's motion and 7„ = 7o(n-l-l) is 
a phenomenologically decoherence rate pO| ]. The peculiar exponent 0.7 in 7„ renders the previous physical explanation 
insufficient. More involved schemes to account for this exponent can already be found in the literature |l^, Here 
we propose a new one based on the previously introduced random evolution schemes. 

The main problem attains the peculiar n-dependency of the argument of the exponential decaying function. In the 
mathematical realm the necessary flexibility comes from a combination of stochastic calculus and the spectral theorem 
for the evolution operator and in the physical one from realizing that not all energy levels of the COM mode can 
be equally affected by a stochastic perturbation. This idea, in a different context in which the trap is coupled to a 
boson reservoir to account for the detected decoherence, has already been paid attention p3[ . 

Let's s tart by considering the spectral decomposition of the evolution operator generated by the Hamiltonian 
(IV. 35) when the laser is tuned to the first blue sideband, i.e. when the Hamiltonian is given by 



Then the evolution operator will be decomposed as follows: 



(IV.41) 



Ur{t) 



OO 

E 

n=0 



(IV.42) 



where e„ = ± rjVl^/n , |e. 



= l/\/2 (| — , n — 1) ± |-|-, n)) {n > 1) and jej) = | + ;0) denote the eigenval ues an d 
eigenstates of ([V.41) respectively and P™^ — |e™^)(e™^| is the projector-valued measure associated to ([V.41). 
The stochastic promotion is performed by substituting e^'t e'!^'t + jQv'^'{s)dB^'{s) in ( IV.42 ) and calculating 



Pint{t) = E f^l"*and(*)p(0)(^/"*and(0)^ • Here uote that the different energy levels are perturbed in a distinct fashion 
determined both by the deterministic functions u™^ (t) and the standard real Brownian motions B^' {t) . The latter 
show correlation properties expressed by the functions g^^m^' (t) '■ 
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dB^'{t)dB'^'{t)=gZ'^^^{t)dt 
The density operator in interaction picture will then be given by 



(IV.43) 



Pint{t) = 



oo 

E ' 

n,m— 
ruz ,m, — ± 



-i Jo («)-/o i'>)dBZ- (s) 



(IV.44) 
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FIG. 2: P-{t) for the orthodox and the stochastic formahsms with p(0) = 
70 = 11.9kHz and 77 = 0.202. See pol. 



, n = 0)( — ,71 = 0| and parameters fi = 470kHz, 



The expectation value in ( IV.44 ) can be calculated with the same techniques as before (cf. also appendix^ and 
drives us to: 



Pi7it 



-i{e'^^ —em" )t — 



(IV.45) 



n,m— 
/ 

rriz ,m, — ± 



where An,^"'^it) = Xn^'^T^ {s)ds with A™^'"-(t) = 
2<^(i)C'-(i)g"™"'^(i). 



The expression ( [V.45 ) already contains the necessary ingredients to arrive at the detected behaviour, since if 
the ion trap is initially set in a Fock state for the COM mode and the ground state for the internal levels, i.e. 
p(0) = | — then the probability p™"''(t) in this scheme is 



j K+i.n+lit) \ , j A„;i.„^ift) 
exp < > + cxp < 



exp <^ -i(e++i - e„+i)t 



exp<^-i(e„+i-e++i)t- 



(IV.46) 
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FIG. 3: P-{t) for the orthodox and the stochastic formahsms with p(0) = | — )(— | (S> pT and parameters = 470kHz, 70 
11.9kHz, r; = 0.202 and Uphotons = 1.5. See Eol. 



Before making physical assumptions let us notice that since the Brownian motions are standard, gn^n^\t) = 1 for 

all n (the brackets mean that both superscripts must be equal) and thus Ai^ri^\t) = for all n. Now we pose the 
most important physical hypothesis, namely the stochastic perturbation depends exclusively upon the energy level of 
the COM mode (at least up to the order of detection we are nowadays capable) . As a first consequence we then can 

claim t hat g^'^it) = 9n'n(^) fo'' ^^'^^ then xl^n'^\t) = and Al^n'^\t) = and also A+'~(i) = A~'+(t), hence 

([V.46) reduces to 



1 + e 2 cos((e++i - e„+i)i) 



K+l r.+ lt*' 

1 + e ^ cos{2r]fltVn + 1) 



(IV.47) 
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FIG. 4: P-{t) for the orthodox and the stochastic formahsms with p(0) — „ <5nQm| — , , m| and parameters Q 
470kHz, 70 = llMHz, 77 = 0.202 and n = 1.5. See ||]. 



Second since for fixed n the internal levels are equally affected, we can also write gn'ni^) — 9n'n{t) ^'^^ each 
n. Finally instead of discussing upon absolute energy values, it is physically more reasonable to talk about energy 
differences and we propose that the stochastic perturbations be introduced in such a way as to have 



e. 



+ -e-^e+-e-+r'^'n% (IV.48) 



where d is an arbitrary exponent, T a constant and where we have assumed for simplicity that the random perturbation 
is a white noise. Note however that it is possible to use more general expressions. Notice the different behaviour of the 
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added term for each distinct subspace of constant COM energy in agreement with the ph ysical hypothesis assumed 



above. Under these hypotheses v^{t) — (t) — F^/^ti'* and after elementary calculations ( IV. 46 ) finally reduces to 



= ^ (l + e^^("+i)"' cos (2ryf7iV^m)) (IV.49) 

This expression shows a clear resemblance to P'^^{t) written above. We believe that both T and d depends 
sensitively upon the particular physical system under study. 

For completeness we also include the expression for p^andj-^-j when the COM mode has an initial state with diagonal 
density- matrix elements P„ : 



prand^^^ = 1 [ 1 + ^ P„e- COS (277^7^^) ) (IV.50) 

This expression has the same structure as the experimental ones shown in [ pO[ . The whole scheme can obviously be 
applied to the carrier, first red sideband and successive excitation too. We include in figs. H, ^ and ^ the predictions 
in the orthodox and the above formalism in the cases when the COM mode is in a Fock state, thermal state and 
coherent state and the internal state is the ground state. 



V. DISCUSSION 



The use of stochastic methods in Hilbert space is of course not new (cf. e.g. ^). The idea of representing open 
quantum systems by means of stochastic processes already appeared in the literature some years ago [^6| and it has 
been widely used in Quantum Optics [ p7| and in the Foundations of Quantum Mechanics Here we pursue the 
line initiated in stepping forward by randomizing not just the (thus stochastic) state vector of the open quantum 
system, but its evolution operator. We find at least three reasons to do that. Firstly when one write a random 
evolution equation for the state vector (thus an Ito stochastic differential equation) an extra term must be included, 
namely the Ito correction. Consider for example the following evolution 



d\-<Pt) = -iH\^t)dt - iA\^jjt)dBt - ^A^\iJt)dt (V.51) 

where the operator A commutes with the hamiltonian H (just for simplicity). From a physical point of view we 
find little intuitive the origin of the Ito correction term, which however appears in a natural way by applying Ito's 
formula to the evolution operator with the stochastic modification Ust{t) = ^-iHt-iABt _ Secondly the use of random 
evolution operators emphasizes the idea that it is the evolution which is random and there is nothing random about 
the Hamil tonia n (and thus the energy levels of the system) , something which may misleadingly be understood from 



equation ( \^.51 ). Finally the use of operators rather than just state vectors opens the possibility of trying to employ 
group-representation techniques [ p9[ and thus of rooting the random nature of the evolution upon possible stochastic 
symmetries. 

This proposal is not intended to solve the so-called macroobjectification problem (better known as the measurement 
problem) by means of a random dynamical reduction process. Indeed it can be readily shown that the state driven 
by a random evolution operator Ust{t) — e^*^*~*^®' is never reduced in clear contrast to these models (see pO|). 
For the case of the previous random evolution operator this can be readily proven. Let Ust{t) — e~*^*~''^^* be the 
random evoluti on op erator of a quantum system ([A, H] = for simplicity). Then the Ito differential equation for the 



state vector is (V.51). Now to check whether this evolution produces dynamical state-reduction or not it is sufficient 



to study the stochastic process (cf. |M) 



where At = Since ^(V'tlV't) = and | (V'olV'o) P = I7 the random evolution is miitary almost surely and then 



dVt = Q^Vt = Vb 



(V.53) 
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Thus there is no actual state-vector reduction process around the eigenvectors of A. This of course differs radically 
for the behaviour of j-f/'t) in stochastic state vector reduction models, where the evolution equation is typically written 
as II, 



dli^t) = -iH\yjt)dt ~liA- Atflijt) + iA~ At)\iJt)dBt (V.54) 

The nonlinear terms appear as a consequence of normalizat ion c onditio ns [E l ] and play no significant role in the 
reduction process. Note the singular difference between eqs. ( V.51 ) and ( V.54| ): despite the fact that both of them 



produce the same master equation (already noted in |42[), only the second one ensures a reduction process taking 
place and this is because of the i factor appearing with the Wiener differential dBt ■ It is an open question whether 
there exists a physical process or not introducing this phase factor in the evolution equation for {ipt}- 
The great utility of stochastic processes to account for the decoherence suffered by a quantum system is its versatility 
to also account for possible intrinsic decohering effects. By this we mean not a fundamental modification of quantum 
principles, as e.g. in |l^, |l^, ^ but just the idea that when a quantum system is described (e.g. an atom interacting 
with the electromagnetic field) some approximations must necessarily be made to be able to analytically handle with 
it and we claim that some of these approximations may hidingly induce decohering effects upon the approximated 
model. In this sense this decoherence can be called intrinsic since no environmental effect is taken into account. 
Undoubtedly this does not deny in any way the possibility of having a system decohered by its environment. Notice 
the special relevance of such a hypothetical effect in quantum systems designed to implement quantum-computational 
and quantum-information-processing tasks, since they usually possess certain degree of complexity which forces us 
to seek for adequate approximations to describe them. A possible relationship of this intrinsic decoherence with 
scalability of quantum-informational and quantum computational systems would also have important consequences 
to find robust mechanisms to process information in a quantum way. 

Besides possible new physical interpretations supporting the use of stochastic processes in Quantum Mechanics, its 
utility to solve certain ME as shown above justifies the search to extend the main result reported here. Currently the 
way to drop the condition of selfadjointness of Lindblad operators is under study. 

Finally a mathematical remark should be made regarding the proof of the previous result. One may wonder why the 
same procedure as in the commuting case, i.e. adding a stochastic term U(t) e~'^*~^ '^o cr(s)dBs ^ jg j^q^ used in the 
noncommuting case. The reason jointly rests upon the Ito's formula and the lack of some derivatives in the case of 
noncommuting operators. Ito's formula |l^ basically states that if the real stochastic process Xt satisfies the Ito 
SDE dXt = a{t,Xt)dt + b{t, Xt)dBt, then the real stochastic process defined as Yt ~ ijj{t,Xt) satisfies the Ito SDE 
given by 



d^it,Xt) = i^ait,Xt)dt^it,Xt) + -b^it,Xt)d^,i,it,Xt)j dt + h{t,Xt)d^i:{t,Xt)dBt (V.55) 

This means that both the first and second partial derivatives of ^{t, x) must exist to be able to apply this formula. 
If the stochastic process is operator-valued as e.g. Q-'^*^-'^^tV ^ then Ito's formula cannot be directly 

applied since the partial derivates of ^-'^^H-ixV cannot be found. This difficulty has been circumvented by changing 
to the Heisenberg picture before introducing the stochastic modifications. However it would be desirable to have an 
Ito's formula for operator-valued stochastic processes valid both for commuting and noncommuting cases. 



VI. CONCLUSIONS 



We have proven how any Lindblad evolution with selfadjoint Lindblad operators can be understood as an aver- 
aged random evolution operator. The proof included here allows us to extend the previous result to nonMarkovian 
situations, though keeping the Lindbladian structure of the master equations, and as a result we have provided a 
straightforward method to solve this kind of master equations. The conjunction of stochastic methods and the spec- 
tral representation theorem has also allowed us to generalize intrinsic decoherence models already present in the 
literature. The mathematical versatility of stochastic methods has also permitted us to propose a generalization to 
the Jaynes-Cummings model, and compare its predictions with experimental results in cavity QED as well as in ion 
traps. We have argued that the main physical advantage stems from the possibility of studying decohering effects 
not necessarily rooted on the environmental action upon the system. A comparison with dynamical collapse models 
reveals that the random unitary operators used here do not produce any kind of state vector reduction. 
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APPENDIX A: SOME USEFUL IDENTITIES 

Though they are elementary we include some useful relations in order to render the text self-contained. First 
we will show how the moments of arbitrary order n of the stochastic process J* v{s)dBs -v{t) a real function- 
are calculated. Let us denote the stochastic process and their nth moments respectively as Xf = v{s)dBs and 
ft 1 " 

/3n{t) = E /p v{s)dBs ■ It is evident that = 0. Xt satisfies the Stochastic Differential Equation dXt — v{t)dBt- 
Applying Ito's formula ||l^ to f{Xt) with f{x) = x" and taking expectation values one readily arrives at 

(2nV 

P2n{t) = ^A"(0 (A.la) 

P2n+l{t)^0 (A.lb) 

where X{t) = f^v'^{s)ds. 

Using the well-known formula cos (A + B) — cos A cos S — sin A sin _B, developing the trigonometric functions into a 
power series and using the previous result one can immediately arrive at 

E[cos{b{t)+ f v{s)dBs)]^e~^ cosb{t) (A.2) 
Jo 



To calculate E[cos^ ZJ, it is convenient to use cos^ A = and then apply (A.2). 

Finally theorem 3 of is quoted with the notation introduced in section 

Theorem 1 If A and B are two fixed noncommuting operators and ^ is a parameter, then 

e«^Be-«^ = e^'^-^ [S] (A.3) 



See original reference for the proof. 
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